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19.  ABSTRACT  (Continued) 


Such  channels  are  known  to  be  rapidly  non-stationary  in  carrier  frequency,  time  or  both. 
This  paper  presents  an  approach  to  modeling  such  doubly  non-stationary  channels.  First 
the  use  of  mutual  coherence  functions  for  the  statistical  representation  of  time-variant 
channel  system  functions  is  briefly  reviewed.  This  is  followed  by  the  presentation  of 
new  time-  and  frequency-separable  forms  of  the  mutual  coherence  functions  that  admit 
doubly  non-stationary  channel  models.  The  WSSUS  model  is  then  shown  to  be  a  degenerate 
case  of  this  more  general  model. 


SUMMARY 


The  standard  approach  to  the  modeling  of  linear  time-variant, 
frequency-selective  communication  cnannels  has  Deen  reviewed  and  extended  to 
allow  for  statistical  variation  in  Doth  time  and  frequency.  One  particular  model 
that  allows  the  channel  statistics  to  vary  in  a  piece-wise  linear  fashion  over 
time  and  frequency  has  Deen  addressed  in  detail.  The  theoretical  aspects  of  the 
model  have  Deen  emphasized.  A  discussion  of  the  use  of  the  DNS  model  for 
practical  communication  system  analysis  can  De  found  in  Reference  2. 
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SECTION  I 
INTRODUCTION 


The  modeling  of  frequency-selective  fading  linear  channels  has  been  a  topic  of 
researcn  for  several  decades,  much  of  it  based  upon  the  work  of  Bello  (Reference 
I).  These  linear  channels  can  be  represented  by  baseband  equivalent 
input-output  system  functions  such  as  a  time-variant  impulse  response  or  a 
time-variant  transfer  function.  The  non-specular  component  of  the  channel 
system  functions  is  usually  viewed  as  the  output  of  a  random  process,  and  the 
random  variations  of  the  channel  can  be  described,  quantitatively  by  the 
statistical  parameters  of  the  associated  random  process. 

When  analyzing  the  effects  of  fading  on  most  communication  systems,  one  can 
assume  that  the  statistical  description  of  the  channel  is  invariant  in  time  and 
frequency.  This  assumption  leads  to  Bello's  so-called  WSSUS  channel  model,  to 
be  discussed  below.  However,  some  communication  systems  operate  over 
propagation  media  with  fading  statistics  that  change  rapidly  in  time  (  e.g. 
meteor  burst  modem)  or  vary  in  frequency  over  the  bandwidth  of  the  transmitted 
signal  (  e.g.  a  direct  sequence  spread-spectrum  HF  modem).  The  WSSUS  modei  is 
inappropriate  for  the  analysis  of  such  systems.  This  paper  presents  a  douDly 
non-stationary  (DNS)  channel  model  to  accommodate  systems  such  as  these,  in 
the  following,  we  first  review  Bello's  general  approach  to  statistical  channel 
modeling,  and  then  derive  time  and  frequency  separable  forms  of  the  channel 
statistics  that  allow  the  channel  model  to  propagate  no  energy  outside  an 
arbitrary  region  in  time-frequency  space.  We  then  show  that  a  weignted 
combination  of  an  arbitrarily  large  number  of  such  models  forms  the  desired  DNS 
model.  Lastly,  we  discuss  one  particular  weighting  that  causes  the  channel 
statistics  to  vary  in  a  piece-wise  linear  fashion  in  time-frequency  space. 


SECTION  2 

GENERAL  BACKGROUND  MATERIAL  ON  CHANNEL  MODELING 


2.1  LOWPASS  EQUIVALENT  CHANNEL  REPRESENTATION. 

The  standard  complex  baseband  signal  representations  will  be  used  to  express 
transmitted  and  received  waveforms.  The  transmitted  signal  is  denoted  z(t)  in 
the  time  domain,  and  2(f)  in  the  frequency  domain.  w(t)  and  w(f)  denote  the 
time  and  frequency  domain  representations  of  the  received  signal.  Both 
transmitted  and  received  signals  are  referenced  to  some  fixed  (angular) 
frequency  oo0  (in  units  of  radians/second),  though  2(f)  and  w(f)  need  not  be 
centered  about  that  frequency.  The  signal  actually  sent  by  the  transmitting 
terminal  can  be  expressed  in  the  time  domain  by  either  side  of  the  following 
equation. 

Real{  z(t)  eJ^ot  }  =  l2(t)cos  oo0t  -  j Qz(t) sino>0t  (l) 

The  complex  waveform  z(t)  equals  I2(t)+  j Qz(t).  where  iz(t)  and  Qz(t)  are  by 
definition  real  waveforms  and  j  denotes  the  square  root  of  -1.  The  signal 
actually  seen  at  the  receiving  terminal  can  similarly  be  expressed  m  tne  time 
domain  as 

ReaK  w(t)  eJa)ot  }  =  Iw(t)cos  ou0t  -  jQw(t) sinco0t  (2) 

where  w(t)  equals  lw(t) ♦  jQw(t),  and  again  lw(t)  and  Qw(t)  are  real  waveforms. 
The  channel  model  generates  the  received  signal  by  passing  the  complex 
transmitted  signal  through  a  time-varying,  linear  system,  as  shown  below. 


Figure  i.  Block  diagram  of  complex  baseband  channel  model. 
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8el!o  in  Reference  '  aefines  tweive  system  equations  that  express  tne  output, 
e't.ner  w(t)  or  w(f).  n  terms  of  v-e  '"put,  e,t’"er  ztf,  or  Z(').  -our  of  tnese  are 
used  m  the  secue.,  anc  are  reproduced  n  tne  'ci.ow  ng  ecuaticns: 
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(4) 


W(f) 
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(6) 


The  variable  t  is  time  in  units  of  seconds,  f  is  frequency  offset  from  the 
nominal  center  frequency  (oo0/2tt)  in  units  of  hertz,  £  is  delay  in  units  of 
seconds,  and  v  is  doppler  frequency  shift  in  units  of  hertz. 

Be l to  named  these  four  system  functions  and  gave  them  physical  interpretations. 
~he  Input  Deiay-Spread  Function.  g(t,£),  is  the  onannei  output  at  time  t  -esuiting 
'rom  the  excitation  of  the  channel  to  a  impulse  £  seconds  oefore  t.  "rhe  Ou.tout 
Doppler-Spread  cunction,  G(f.v).  is  the  spectral  response  v  hertz  above  f, 
resulting  from  the  excitation  of  the  channel  to  a  frequency  impulse  (i.e.,  time 
domain  cissoid)  at  f.  The  Time-Variant  Transfer  Function,  T(f,t),  is  the  transfer 
function  m  the  frequency  variable  f  at  time  t.  Finally,  suppose  the  transmitted 
signal  is  decomposed  into  infinitesimal  elements  that  are  first  delayed,  and  then 
doppler  shifted  before  being  linearly  combined  to  form  tne  channel  output.  Then 
the  Delay-Doppler  Spread  Function,  u(£,v),  is  the  complex  weignting  given  to  the 
component  delayed  by  £  seconds  and  shifted  by  v  hertz. 

Reference  l  shows  that  the  four  system  functions  are  related  by  Fourier 
transforms  in  one  or  both  arguments.  The  various  reiationsmps  are 
mathematically  expressed  by  Equations  7  through  18  and  illustrated  by  Figure  2 
(where  arrows  point  in  the  direction  of  the  Fourier  transforms)  . 


g(t,£) 
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gCt.O 

=  |  J  G(f.v)  e  J2ft(vt<f)  dv  df 
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T(f.t) 

=  J  g(t.O  e  ~ j2ttCt 
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T(f.t) 

=  |  |  Utt.v)  e  J27t(vt-Cf)  dv  dC 

(12) 

U5.V) 

=  g(t.O  e  ~ j2rcvt  dt 

(13) 

U(£.v) 

=  [  G(f.v)  e  J2rc^f  df 

(14) 

U(C.V) 

=  J  |  T(f.t)  e  j2Tt(^f-vt)  df  dt 

(15) 

G(f.v) 

=  J  T(f.t)  e  *j2«vt  dt 

(16) 

G(f.v) 

=  J  U(£,v>)  e  -J2^f  dC 

(17) 

G(f.v) 

=  |  \  g(U)  e  -j2*(vf£f)  <jt  d C 

(18) 

2.2  STATISTICAL  FORMULATION  OF  RANDOM  CHANNELS. 

Random  channels  can  be  characterized  by  statistical  descriptions  of  any  one  of 
these  four  system  functions,  it  is  common  to  assume  the  real  and  imaginary 
components  of  the  system  functions  are  independent  and  have  Gaussian, 
zero-mean  first  order  statistics.  The  stochastic  properties  of  the  random 
system  functions  are  then  completely  characterized  with  the  addition  of  the 
second  order  statistics.  The  second  order  statistics  are  expressed  by  the 
correlation  properties  of  the  system  functions  in  the  form  of  what  we  term 
“mutual  coherence  functions".  The  mutual  coherence  functions  of  the  four 
system  functions  are  defined  as  follows. 


RgCt.SlE.T)) 

=  g(t.E)M  gCs.-q) 

09) 

Rj(f ,  1 :  t.  s) 

=  T(f.t)"  T(l.s) 

(20) 

.  -  .  U  ' 

Ro(f,  liv.ji) 

=  G(f.v)w  G(l.ji) 

(21) 

a M&ajui 

39 rw' 

BMM 

RuCE.^iv.m) 

=  U(£.vf  U(rt.M) 

(22) 

^  IJ  ^  J  W  J 

where  *  denotes  the  complex  conjugate  operation  and  the  overbar  indicates  the 
ensemble  averaging  operation.  As  was  the  case  with  the  system  functions,  from 
any  one  of  the  mutual  coherence  functions  one  can  derive  the  other  three.  They 
are  interrelated  by  double  or  quadruple  Fourier  transforms  as  discussed  in 
References  1,  expressed  mathematically  by  Equations  23  through  34  and  shown 
pictorial ly  in  Figure  3. 
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Figure  2.  Fourier  transform  relationships  between  the  system 
functions  g  ,  T  .  G  .  and  u  . 
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Rg(t.s;£.Ti)  =  I  j  RT(f,l;t.s)  e  j27TCn !-^r)  df  dl  (23: 

Rg(t.s;£,Ti)  =  Rf/S.riiv.pO  e  j2ir(jis-vt)  gv  g^  (24; 

Rg(t.s;^,T|)  =  RgCf.liv.ji)  e  j27T(jis-vt+T[t-Sf)  dv  d^  df  g;  (25; 

RT(M:t,s)  =  f  [  RgCt.siS.'n)  e  j2rr(Cf-r(  I)  d£  d<n  (26) 

RjCf.ht.s)  =  R0(f . l;v. jji)  e  j2rr(p.s-vt)  dv  g^  (27; 

RyCf. t:t.s)  =  R|j(^.tl;v,|i)  e  )27T(^s-vt+£f-Ti  l)  dv  d^  d£  dT^  (28) 

J  J  J  J 

R(j(^.il;v,ji)  =  Rg(t,si£,T|)  e  j2rr(vt-|is)  gt  ds  (29) 

RU(C.,r\;V.Jl)  =  f  Rg(f , l;v\ ji)  e  j27T(T[  |-£f )  gf  gi  (30) 

R(j(£,1];V,;i)  =  j  ||  Rr(f,l;t,s)  e  j2TT(T\  ! -£( f ♦  Vt~ JJ.S)  gf  g)  gf  gs  (3  1 ) 

RG(f . =  |  RT(f.l:tTs)  e  j2rr(vt-ps)  dt  as  (32/ 

R<3(f,l;v, ;i)  =  |  R|j(C,T|;v,|i)  e  d-q  (33) 

R(3(f,l;v,jj)  =  Rg(t,s;^,Ti)  e  j27T(vt-jas+£f-T| t)  dt  ds  d£  dT|  (34) 


2.3  WSSUS  CHANNEL  MOOEL. 

The  WSSUS  channel  model  is  obtained  from  the  above  general  statistical  form  by 
assuming  the  first  and  second  order  statistics  are  constant  over  (or  independent 
of)  time  and  carrier  frequency.  The  time  independence  is  called  the  wtde-sense- 
stationary  (WSS)  assumption  and  the  carrier  frequency  independence  is  called 
the  uncorrelated  scatterer  (US)  assumption.  The  effect  of  each  of  the  WSS  and 
US  assumptions  upon  the  form  of  mutual  coherence  functions  is  first  treated 


separately  in  the  discussion  to  follow.  Then  taken  together  the  W5SUS  channel 
follows  immediately. 
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Figure  3.  Fourier  transform  relationships  between  the  mutual 
coherence  functions  of  the  system  functions  g  .  T  .  G 
and  U  . 


The  WSS  assumption  causes  a  simplification  of  the  functional  form  of  all  four 
mutual  coherence  functions  discussed  above.  The  WSS  forms  of  the  mutual 
coherence  functions,  denoted  with  a  tilde  throughout  the  sequel,  are  as  follows. 


Rg(t,s:£.Tl)  =  ffg(s-t;£.T() 

RT(f  ,l;t,s)  =  ST(f.l;s-t) 

RgCf.liv.jj)  =  pQif.Uv)  S(p-v) 

Ry(^.T\:V.n)  =  PyCS.niV)  S(p~v) 


Substituting  s+v  for  t  allows  the  WSS  forms  of  Equations  23  through  34  to  be 
rewritten  as  follows. 


Rg(r;C.i\) 

Rg(r;£:.Ti) 


fiT(f.l;r)  e  j27T(Tll-^f)  df  d| 
PuiK.r[-.v)  e  )27Tvr  dv 

f  pe(f.l;v)  e  j2TT(vr^Kf)  dv  df  d, 


(39) 

(40) 

(41) 


ST(f.l;r) 

PT(f.l:r) 

PT(f.':r) 


Pg(r:^.T\)  e  j2tr(£f-Tii)  d(r  dT^ 
p6(f,l;v)  e  J2™:  dv 
[  PyU^iv)  e  j27t(vr+£f-T\l)  dv  d£  6t\ 


(42) 

(43) 

(44) 


PjZ,r\w) 

Pytt.'HiV) 


Pg(r;^.T|)  e  "j2ttvr  dr 
PG(f,Uv)  e  j2*eqKf)  at  di 
PT(f,l:r)  e  j2tt(ti  Kf-vr)  df  d!  dr 


(45) 

(46) 

(47) 


?GCf-1:v;) 


PG(f.l;v) 


PT(f,l;r)  e  ~J2TTvr  dr 
=  |  Py(^,T[;v)  e  j2Tt(£f-Til)  d£  dTj 

=  |||  Pg^C.Ti)  e  j2ft(-vr^f-Til)  dr  d£  dm 


(48) 

(49) 

(50) 


The  above  equations  show  that  the  mathematical  relationships  between  the 
non-singular  factors  of  the  WSS  forms  of  the  mutual  coherence  functions,  like 
the  general  forms  themselves,  are  related  by  Fourier  transforms.  Figure  4 
shows  these  Fourier  transform  relationships.  The  non-singular  factors  in  the 
WSS  forms  of  the  mutual  coherence  functions  have  the  following  physical 
interpretations. 
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ftg(r;£.Tj)  Cross  correlation  of  pair  {g(t,C).g(t.ri )]  for  fixed  £  and  -q, 

3y(f.l:r)  Cross  correlation  of  pair  [TCf .t).T(l.t)]  for  fixed  f  and  I. 

p0(f .  I ; v»)  Cross  spectrum  of  pair  tT(f.t).Td.t)]  for  fixed  f  and  I. 

PgCf.  f  ;v)  as  a  function  of  v  is  the  power  spectral  density 
of  the  fluctuations  of  T(f.t)  in  the  variable  t. 

Cross  spectrum  of  pair  [g(t.£).g(t.^)]  for  fixed  £  and  r\. 
P\j(£,£:v)  as  a  function  of  v  is  the  power  spectral  density 
of  the  fluctuations  of  g(t,0  in  the  variable  t. 

The  WSS  assumption  affects  the  nature  of  the  system  functions.  G(f.v)  and 
U(£,vO  have  the  character  of  non-stationary  white  noise  in  the  variable  v.  g(t,£) 
and  T(f,t)  have  the  character  of  stationary  colored  noise  in  the  variable  t. 

The  US  assumption  causes  a  different  simplification  in  the  forms  of  the  mutual 
coherence  functions.  These  US  forms,  denoted  with  a  circumflex  throughout  the 
sequel,  are  as  follows. 


Rg(t.s;£.-q)  =  Pg(t.s;^)  S(^-rl) 

(51) 

nT(f.l:t.s)  =  ftT(l-f;t.s) 

(52) 

R  q(  f,  l :  V,  jji)  =  ftG(l-f;V.>l) 

(53) 

RutS.rpv.ji)  =  P>u(S;v.^i)  8(£-q) 

(54) 

Substituting  Q  for  l-f  allows  the  US  of  Equations  23  through  34 
as  follows. 

to  be  rewritten 

Pg(t,s;£)  =  [  AT(Q;t,s)  e  J2Tr^Q  dQ 

(55) 

Pg(t,S;E)  =  e  j27T(j-lS-Vt)  dv  d^ 

(56) 

Pg(t.s;0  =  |  &g(Q:V-M)  e  j2rr(MS-vf^Q)  dv  d^  dQ 

(57) 
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Pu(£.t];v)  S(m-v) 


S T ( r .  i :r) 


PgCt.s:0 

Pg(t,S;0 

pqtt.s£) 

ftjCQit.s) 

ftT(Q:t.s) 


PG(f,  I :  v)  S(|i-v) 

Fourier  transform  relationships  Detween  the  wSS 
'arms  of  the  mutual  coherence  'unctions. 


ftT(Q:t.s)  e  J27TCQ  dQ 
j  P\Jtt-.v,n)  e  )2Tt(MS-vt)  dv  dM 
J  Pg(Q;v,^l)  e  j2TT(|is-vt*£Q)  dv  dju  dQ 


=  |  Pg(t.s;0  e  -i2trCQ  d£ 


AG(Q;v,|i)  e  ^(us-Vt)  dv  d^ 


e  j2TT(Ms-vt-£Q)  dv  dp  d(r 


"/ ^  "r  -  ■' 
•v/A-  *. 


i\j(£;v.*i) 

=  j  £g(t,s;£)  e  j2TT(vt-|is)  dt  ds 

(60 

PuUuv.m) 

=  f  AG(Q;v.^)  e  l27t^  dQ 

(62) 

=  Aj(Q:t.s)  e  j2fT(£Q*vt-Ms)  dQ  dt  ds 

(63) 

A<3(Q:v.ji) 

=  |  [  ftT(Q:t.s)  e  j 2tt( vt — jis)  dt  ds 

(64) 

R<3(Q;V,ji) 

=  l\,(£;v.ji)  e  J2tt-^Q  d£ 

(65) 

Aq(Q;V.|i) 

=  [  J  J  Pg(t.s;£)  e  j27t(vt-]is-£Q)  dt  ds  d£ 

(66) 

The  Fourier  transform  relationships  between  the  US  forms  of  the  various  mutual 
conerence  functions  are  shown  Dy  Figure  5.  The  non-singular  factors  in  the  US 
forms  of  the  mutual  coherence  functions  have  the  following  physical 
interpretations. 

?g(t.s;£)  Cross  spectrum  of  pair  [T(f,t),T(f,s)]  for  fixed  t  ana  s. 

3g(t.  t;£)  as  a  function  of  £  is  the  power  spectral  density  of 
the  fluctuations  of  T(f.t)  in  the  variable  f. 

R-r(Q;t,s)  Cross  correlation  of  pair  [T(f,t),T(f,s)l  for  fixed  t  and  s. 

Cross  correlation  of  pair  [G(f,v),G(f,ji)]  for  fixed  v  and  ji. 

Cross  spectrum  of  pair  [G(f.v),G(f.ju)]  for  fixed  v  and  jjl 
PyCCiv.v)  as  a  function  of  £  is  the  power  spectral  density  of 
the  fluctuations  of  G(f,v)  in  the  variable  f. 

The  US  assumption  affects  the  nature  of  the  system  functions.  g(t,£)  and  U(£,v) 
have  the  character  of  non-stationary  white  noise  in  the  variable  £.  G(f.v)  and 
T(f.t)  have  the  character  of  stationary  colored  noise  in  the  variable  f. 

The  WSS  and  US  assumptions  taken  together  is  called  the  WSSUS  assumption. 
The  WSSUS  form  of  the  mutual  coherence  functions,  denoted  with  an  over-bar 
through-out  the  sequel,  are  given  below. 
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Figure  5. 
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Fourier  transform  relationships  between  the  US 
forms  of  the  mutual  coherence  functions. 


(67) 

(68) 

(69) 

(70) 
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Substituting  |*Q  for  f  and  s+r  for  t  allows  the  WSSUS  forms  of  Equations  23 


tnrougn  34 

to  be  written  as  follows. 

Pg(t;£) 

=  [  Rj(Q:r)  e  j2tt^Q  dQ 

(7i; 

Pg(r;E) 

=  j  Pu(c;v)  e  J2"hvr  dv 

(72; 

P n(r:0 

=  f  f  R"g(Q;v)  e  j2Tt(vr*£Q)  dv  dQ 

(73; 

Rj(Q;r) 

=  j  Pg(T;S)  e  -J2tt^Q  d£ 

(74) 

=  PqCQiv)  e  J27tvr  dv 

J 

(75) 

=  [  |  PyCSiv)  e  j2Tt(vr-^Q)  dv  d£; 

(76) 

P[j(S;v) 

=  Pg(r;C)  e  ")2Ttvr  dr 

(77) 

Pi[/£:v) 

=  PG(Q;v)  e  ]2tt£;q  dQ 

(78) 

Pu(^;v) 

=  fty(Q:T)  e  j2Tt(£;Q-vr)  dQ  dr 

(79) 

ftgCQiv.y) 

=  Py(Q;t,s)  e  j2rc(vt-ps)  d^  ds 

(80) 

Aq(Q:v.ji) 

=  f  PuCCiv.ji)  6  J2tt-CQ 

j 

(81) 

PG(Q;  v\jl) 

-  j  Pg(t,s;C)  e  j2Tt(pt-ys-cQ)  gt  ds  d£ 

(82) 

The  Fourier  transform  relationships  between  these  are  shown  schematically  in 
Figure  6.  The  non-singular  factors  in  the  WSSUS  forms  of  the  mutual  coherence 
functions  have  the  following  physical  interpretations. 

Pg(r;£)  Delay  cross-power  spectral  density.  P'g(Oi^)  is  the 
Delay  power  density  spectrum. 

Ry(Q;r)  Time-frequency  correlation  function. 

PqCQiv)  Doppler  cross-power  spectral  density.  Pq(0;v)  is 
the  Doppler  power  density  spectrum. 

P(j(£;v)  Scattering  function.  The  power  spectral  density  in 
delay  and  Doppler  of  T(f,t) 


The  W55US  assumption  affects  the  nature  of  the  system  functions.  g(t,£)  has 
the  cnaracter  of  non-stationary  white  noise  in  the  variable  £  and  stationary 
colored  noise  in  the  variaDle  t.  U(£,v)  has  the  cnaracter  of  non-stationary  white 
noise  in  both  variables.  G(f.v)  has  the  cnaracter  of  non-stationary  wmte  noise 
in  the  variable  v  and  stationary  colored  noise  processes  in  the  variable  f.  T(f.t) 
has  the  character  of  stationary  colored  noise  in  both  variables. 


;  v)  S(£-g)  S(pi-v)  R  T  (Q  ;  r ) 

I  I 

I  1 

\  / 


P^(Q;v)  o(u-v) 


Figure  6.  Fourier  transform  relationships  between  the  WSSUS 
forms  of  the  mutual  coherence  functions. 


SECTION  3 

SEPARABLE  FORMS  OF  THE  MUTUAL  COHERENCE  FUNCTIONS 


3.1  TIME-SEPARABLE  FORMS. 

Recall  from  Equations  37  and  38  tnat  the  WSS  assumption  leads  to  a  singularity 
of  the  form  S(jjl-v)  in  the  mutual  coherence  functions  R0  and  R^.  Here  we 
consider  replacing  8(yrv)  in  Equations  37  and  38  with  seme  arbitrary 
non-singular  factor,  denoted  P(|i-v).  Then  we  snow  tne  corresponding  effects 
upon  Rj  and  Rg.  We  start  with  a  pair  of  functions,  say  £(■)  and  f(-).  tnat  are 
related  Dy  the  following  one-dimensional  Fourier  transformation. 

P(x)  =  |  S(s)  e  -J2ttxs  ds  (83) 

From  Equation  27  we  have 

R-(f,;;t,s)  =  |  f  Rrjvf.hv.jj)  9  :2"r(us-vt)  dv  du 

-resuming  R^f.iiv,^)  is  :n  the  form  we  can  wr-te 

r  r  r  1 

Rj(f,l;t,s)  =  J  pQ(f,l:v)  j  e  J2tt( j_is- vt)  j  dv 

L  -i 

Replace  jus-vt  with  v(s-t)  *  (pi-v)s  to  get  the  following  equation. 

Rj(f,l;t.s)  =  Pe(f,l;v)  e  j2ftv(s-t)  f(jj-v)  e  j2TtCi_i.- v)s  du  dv 
Substituting  v+f3  for  p  and  dfi  for  dyi  yields 

RT(f.l;t,s)  =  j  PG(f,l:v)  e  j2rcv(s-t)  dv  j  e  j2tt3s  d£  (87) 

From  Equation  43  and  63  we  get 

Rr(f.i;t.s)  =  3T(f,l;s-t)  £(s)  .  (88) 


Therefore  replacing  the  singularity  S(p-v)  in  the  WSS  form  of  Pg  with  some 


'’j-  r  'ir 


non-singular  function  F(^i-v)  introduces  the  solely  time-dependent  factor  2f(s) 
into  tne  WSS  'orm  of  R;.  Equation  88  could  be  rewritten  as 


R~(f,!;S-tr,s) 

=  ST(f.l;r)  ?(s)  . 

(89) 

so  that  all  of  the  s,  or  time,  dependence  has  been  separated  out  of  St  and  is  in 
?(s).  a  similar  treatment  relates  tne  time-separable  ana  WSS  forms  of  the 
'emaining  mutual  conerence  functions.  The  following  equations  snow  all  of  tne 
t  me-separaoie  'orms  ana  should  be  compared  directly  with  Equations  35  tnrougn 
35. 

:v  ■  i  * 

=  Sg(s-t;C.m,)  2f(s) 

(SO) 

^  ,  .  .  -.0  / 

=  ST(r.:;s-t)  ?(s) 

(9  0 

=  3'/.r.v,y) 

=  S3(f.l;v)  f(jj-v) 

(92) 

-  • ! ,  ■"  .  v .  __ 

~~e  vSS  ■'  sr~'  t 

"er  s  ust  a  special  case  ;f  the  t 

'r^e-seo3r3D:e  ^ 

part  cular,  wnen  §(s)  2l  and  {“(ju-v)  =  S(pi-v).  Figure  7  •$  tne  mere  genera, 
t  "'e-separaoie  version  of  Figure  4.  The  functions  Sg  ,  St-  ,  S3  and  r 3  are  stnl 
'elated  mathematically  by  Equations  39  tnrougn  50. 

3.2  FREQUENCY-SEPARABLE  FORMS. 

Fecal i  from  Equations  51  and  54  that  the  US  assumption  leads  to  a  singularity  of 
the  form  S(£-n)  in  the  mutual  coherence  functions  F-  and  Fij.  uere  we  consider 
'epiacmg  in  Equations  15  and  18  witn  some  arbitrary  ncn-s.nguiar 

'actor,  denoted  f(£-ri).  Then  we  show  the  corresponding  effects  upon  Rt  and 
30.  we  start  with  a  pair  of  functions,  say  &(■)  and  f(-).  that  are  related  by  the 
following  cne-dimensional  Fourier  transformation. 


=  |  21(1)  e  ~ J 2rtx l  d| 


fQj) 
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3g(r;£.'n)  §(s) 


PG(f ,  I ;  v) 


cigure  7.  Fourier  transform  "elat  cosmos  oetweer  •me 
time-seoaraoie  ‘'or^s  of  'me  mutual  oorererc9 
functions. 


From  Equation  25  we  nave 


PT(f.l;t,s) 


j  Mg 


Ra(t,s;^.r[)  9  J 2rr(c -m  G£  G 


Presuming  Rg(t,s;C.T[ )  is  in  the  form  Pg(t.s;£)  we  can  //n 


Rr(f .  Iit.s) 


f  Pa(t,s;0  f  |  r(£-m)  e  c- 


L  J 


Peoiace  *f-n!  witn  ^ - m ) !  to  get  the  'cilcwnc  ecuafon. 


Rr(M:t.s)  =  I  P»Q(t.s;0  e -j2^0-f)  \  f7£-m )  e  J^TT^-'n  )l 

S.cstituting  '-3  'or  n  and  -dS  “'or  dm  yields 

=  i  e„(t.s;C)e-J2JrC(l-f)a£  :  f'3 )  e  i2^5 :  03 


=  mM;t,s) 


m 


Pa(t.s:£)  fXt-T\) 

y  I  I  ! 


ftT(Q:t.s)  *(l) 


ftG(Q;v,]i)  5f(l) 

Figure  8.  Fourier  transform  relationships  between  the 

frequency-separaDie  forms  of  the  mutual  coherence 
functions. 

From  Equations  58  and  94  we  get 

RT(f.l;t,s)  =  ftT(l-f:t.s)  £(1)  (99) 

Therefore  replacing  the  singularity  S(£-t[)  in  the  US  form  of  Pg  with  some 
non-singular  function  f(£-T\)  introduces  the  solely  frequency-dependent  factor 
#0)  into  the  US  form  of  Rj.  Equation  99  couid  be  rewritten  as 

RT(1-Q.!;t,s)  =  ftT(Q:t.s)  Sf(l)  .  (100) 

so  that  all  of  the  I,  or  frequency,  dependence  has  been  separated  out  of  Ry  and 
is  in  <if(l).  A  similar  treatment  relates  the  frequency-separable  and  US  forms  of 
the  remaining  mutual  coherence  functions.  The  following  equations  show  all  of 
the  frequency-separable  forms  and  should  be  compared  directly  with  Equations 
51  through  54. 


Rg(t.s;£/q)  =  Pg(t.s;S)  ftS-n) 


R  j(f .  1  :t,s) 


=  ftT(|-f;t.s)  2(1) 


(101) 


(102) 


RG(f,l;v,ji) 


U03) 


=  RGU-f;v.M)  2fCU 

Ry^.TlJV.M)  =  ^(SiV.JA)  (104) 

The. US  form  then  is  just  a  special  case  of  the  frequency-separable  form,  in 
particular,  when  2(1)  *1  and  f(E-T\)  =  S(£-t|).  Figure  8  is  the  more  general 
frequency-separable  version  of  Figure  5.  The  functions  .  ftr  ,  $g  30(3  3rs 
still  related  mathematically  Dy  Equations  25  through  66. 


3.3  TIME-  AND  FREQUENCY-SEPARABLE  FORMS. 


The  derivations  of  the  two  previous  subsections  can  easily  be  drawn  together  to 
provide  the  time-  and  frequency-separable  forms  of  tne  mutual  coherence 
functions  in  terms  of  the  WSSUS  forms  and  the  arbitrary  Fourier  transform 
functional  pairs  [?(•).?(•)]  and  (2(-),f(-)l  .  These  forms  are  snown  in  the 
following  four  equations  and  should  be  compared  with  Equations  67  through  70. 


C)  2(s) 


i  (c-n 


Figure  9. 


Pg(Q;v)  2(l) 


Fourier  transform  relationships  between  the  time-  and 
frequency-separable  forms  of  the  mutual  coherence  functions. 
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RgU.S;£,Tl)  =  Pg(S-t;£)  r(£-T\)  2f(S) 

RT(f,l;t.s)  =  Rj(l-f;s-t)  tf(s)  #(!) 

Rtfr. Isv.ji)  =  FG(|-f:v)  V(^x-v)  *(l) 

r^c.tv.v.^)  =  pu(k-.v)  f(M-v)  rcc-ri) 


u  05; 


(106) 


(107) 


008) 


^he  WSSUS  form  is  just  a  special  case  of  the  time-  and  frequency-separaDle 
form.  Figure  9  is  the  more  general  version  of  Figure  6.  The  functions  Pg  ,  R"y  , 
°0  and  Py  are  still  related  mathematically  Dy  Equations  71  through  82. 


3.4  UTILITY  OF  THE  SEPARABLE  FORMS  FOR  THE  DNS  CHANNEL  MODEL. 

The  general  forms  of  the  four  mutual  coherence  functions,  defined  previously  Py 
Equations  19  through  22,  provide  complete  four  dimensional  statistical 
characterizations  of  the  channel.  These  forms  are  clearly  non-stationary  in 
time  and  frequency.  Put  do  not  lend  tnemseives  to  an  intuitively  appealing,  or 
even  manageaoie.  parameterization.  The  assumption  of  separaoility  has  provided 
the  more  tractaoie.  Put  less  general,  forms  of  the  mutual  coherence  functions 
given  py  Equations  105  througn  108.  These  forms  are  too  restrictive  to  Pe  used 
for  direct  characterization  of  arbitrary  variation  of  the  channel  statistics  in 
time  and  frequency.  However,  these  do  provide  a  formalism  to  turn  the  channel 
off  over  any  region  in  the  frequency-time  space.  This  capability  provides  the 
oasis  for  very  general  non-stationary  channel  models. 

Consider  the  separable  form  of  Ry  given  by  Equation  106.  with  two  variable 
substitutions,  Q  =  I  -  f  and  r  =  s  -  t  .  this  equation  can  be  rewritten  as 

Ryd-Q.lis-r.s)  =  R'-KQir)  ?(s)  £(l)  (109) 

R’yCQir)  specifies  a  WSSUS  channel  model,  and  thereby  contains  the  requisite 
information  to  model  the  channel  over  any  range  of  time  and  carrier  frequency 
that  the  channel  remains  approximately  stationary.  The  ?(s)  &(1)  factors 
provide  the  flexibility  to  drive  the  right  side  of  Equation  109  to  zero  at 
arbitrary  regions  of  time  and  carrier  frequency,  even  when  Q  and  v  equal  0.  A 
value  of  zero  corresponds  to  no  energy  being  propagated. 


ssfiSS 
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In  the  next  section,  an  array  of  channel  models  will  be  combined  to  form  a 
composite  channel.  Each  will  have  distinct  ?(s)  £(l)  factors  to  impede  energy 
propagation  outside  overlapping  regions  of  time  and  carrier  frequency.  The 
FT-rCQir)  factor  for  each  component  model  will  provide  a  WSSUS  approximation  of 
the  channel  statistics  within  the  associated  time-frequency  regions  of  non-zero 
S(s)  ft(l).  This  non-stationary  composite  model  provides  an  intuitively  appealing 
means  to  parameterize  the  channel  --  by  specifying  a  two  dimensional  WSSUS 
characterization,  In  the  form  of  Wy(Q:r),  from  region  to  region  over  arbitrarily 
large  intervals  of  time  and  carrier  frequency. 
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SECTION  4 

THE  DOUBLY  NON-STATIONARY  (DNS)  CHANNEL  MODEL 


Recall  from  Section  2.3  tnat  a  WSSUS  channel  model  can  De  specified  Dy  any  one 
of  the  four  two-dimensional  mutual  coherence  functions.  The  most  intuitively 
appealing  of  these  parameterizations  is  the  scattering  function,  denoted  pJj  in 
Section  2.3,  a  two-dimensional  power  spectrum  in  the  delay-doppler  domain, 
exniDitmg  no  variation  in  time  or  carrier  frequency..  This  section  presents  a 
non-stationary  channel  model  that  allows  gradual  variation  in  the  scattering 
function  with  time  and  frequency.  The  “gradual"  restriction  on  the  variation  of 
the  scattering  function  with  frequency  and  time  allows  the  channel 
representation  over  any  small  region  of  time  and  frequency  to  De  well 
approximated  with  the  WSSUS  formalism.  A  channel  having  this  locally 
stationary  characteristic  is  sometimes  called  a  quasi-wssus  or  QWSSUS 
channel. 

In  this  development  of  the  non-stationary  channel  representation,  direct  use  of 
the  scattering  ^unction  is  inconvenient,  and  we  will  instead  use  its  dcuDle 
Fourier  transform,  denoted  previous  Dy  fify.  In  addition  the  frequency-time 
aomam  is  allowed  to  De  arDitrarily  iarge,  Dut  is  partitioned  into  a  rectangular 
grid  with  the  frequency  granulation  indexed  Dy  m  and  the  time  granulation 
indexed  Dy  n  as  shown  in  Figure  10. 


4.1  COMBINING  MULTIPLE  WSSUS  SYSTEM  FUNCTIONS  TO  OBTAIN  DNS 
SYSTEM  FUNCTIONS. 

we  assume  that  an  independent  random  WSSUS  system  function,  say  the 
time-variant  transfer  function,  is  available  for  each  grid  line  crossing  point. 
Each  such  transfer  function  exhibits  distinct  statistics  completely  specified  Du 
one  of  the  mutual  conerence  functions,  for  example  Wy.  For  this  development, 
we  augment  the  notational  conventions  for  system  functions  and  mutual 
coherence  functions  with  the  grid  point  indices.  Thus  the  time-variant  transfer 
function  at  the  (m,n)3il  grid  point  is  denoted  as  Tn  m(f,t),  and  the  corresponding 
WSSUS  mutual  conerence  function  is  denoted  as  3y  i  n_m(Q,t),  which  in  accord 
with  Equations  8  and  20,  and  the  augmented  notation  tor  T,  is  written  as 
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Rt  |  n,m^*r)  =  Tn  ^Cl.s)  .  (i  10) 

The  time  and  frequency  origin  of  Tn  m(f,t)  is  taken  relative  to  the  associated 
grid  point  coordinates  Tn  and  Fm  .  we  desire  to  combine  the  independent 
time-variant  transfer  functions,  variously  offset  in  time  and  frequency  to  tneir 
associated  grid  points,  in  such  a  way  that  the  resulting  composite  time-varying 
transfer  functions  satisfy  the  following  objectives: 

1.  The  composr  ystem  function  should  be  QWSSUS  over  time-frequency 
regions  tha*  .  small  in  comparison  to  the  local  grid  cell  size. 

2.  In  the  vie*-'  ty  of  a  grid  point,  the  composite  system  function  should 
exhibit  Q  jSUS  statistics  that  match  the  WSSUS  mutual  coherence 
function  corresponding  to  that  grid  point. 

3.  The  QWSSUS  statistics  of  the  composite  system  function  should  vary 
smoothly  between  the  grid  lines. 

4.  if  identical  WSSUS  mutual  coherence  functions  are  specified  for  every 
grid  point,  the  composite  system  functions  should  degenerate  back  to  a 
WSSUS  channel  characterization  spanning  the  entire  grid. 

Tne  composite  system  functions  require  a  few  more  extensions  to  our  notaticnai 
conventions,  we  denote  the  absolute  frequency  vanaDie  as  f3  .  wmen  s  the 
abscissa  of  riqure  50.  and  the  absolute  time  variable  as  t3  .  wmen  ;s  the 
ordinate  of  Figure  10.  These  are  related  to  the  the  grid  point  relative  frequency 
and  time  variables,  t  and  f,  by  offsets  to  the  associated  grid  point  co-orainates 
oy  f  =  fa  -  -m  and  t  =  ta-Tn  .  Also,  system  functions  that  are  not  assumed  to 
nave  WSSUS  statistics  are  expressed  in  bold  type.  This  will  allow,  for  example, 
the  system  functions  Tn,rmr^)  anc!  Tn,m^r-^  -°  3e  defined  differently  oeicw. 

The  above  objectives  can  be  well  satisfied  by  a  composite  time-variant  transfer 
function  of  the  form 

M  N 

T(ra.ta;  =  Z!  wn,m^a‘rm  •  ^a"  1  n>  Tn,m^a'rm  ^a'  n-1 

m-0  n-0 

where  wn  m(fa-Fm , ta-Tn)  is  assumed  to  be  everywhere  real  and  positive,  we 
now  explore  the  suitability  of  this  definition  of  T(f3,ta).  First  defme  the 
weighted,  and  therefore  non-wSSUS.  transfer  function  Tn  rn(f,t}  as 

Tn.m^-t)  =  wn,m(f- *■)  •  (!I2) 


23 


Its  mutual  ccnerence  function  is  tnen  given  Dy 
R  T  I  n,m^'  =  Tnrn(f.t)  TnmU.s)  .or 


=>T  j  n_mU-Q,l;s-r,s)  =  Tn  m(l-Q.s-r'>"  Tn  m(l.s) 


=  wn  rri ( l-Q, s-r)  wnifin(l.s)  R t  |  n,m(^'J) 


"en  Equation  ;  i i  can  De  rewritten  m  terms  of  Tn  m  as  follows. 


in 

T(f3,tg.)  =  S  2  Tn.m^a’^m  •  ta-T"n) 

rn- 0  rK) 


its  mutual  conerence  function  is  given  t)y  the  following  equation. 


(M3i 


(1  Mi 


(116' 


n  N  M  N 


u-wh**  -  Z  I  I  S  To. 

*n«0  ti  w)  n  -0 


m^a'^m^a'  '  0)  T-v  tt'  l^~pm  '•  S3-  rv  1 


n ',.m '  'a  m  '•  3a 


S' nee  tne  time-variant  transfer  functions  associated  with  distinct  grid  points 
are  ndepercent,  the  above  equation  can  be  re-written  as 


■ 

Tf3.l3.t3.S3/  =  y  y  I  n.m(fa'Fm'la"rm;  ^‘‘n-Sa' 

mn)  n-0  1 


C  :s; 


or.  jsmg  the  suDstitution  of  Equation  111  and  replacing  fa  and  ta  witn  la-Q  and 
s3-r,  respectively,  one  arrives  Equation  1  19  Delow. 

9  ytu-Q. l3,sa-t,s3)  = 

4  (119) 

M  N  __ 

S  S  wn.m^a'c:m Sa'^n*7^  wn,m^a-,’m'  Sa'^n)  Rl|n,m^-r^ 

t«0  n«0 

Now  the  character  of  w  needed  to  satisfy  our  oDjectives  can  De  ascertained.  For 
ar.  channels  of  interest,  it  can  be  assumed  that  there  exists  some  Qmax  and  rmax 
called  respectively  tne  maximum  correlation  frequency  and  the  maximum 


correlation  time,  for  whicn  Ry  |  nim(Q.r)  is  essentially  zero  when  |  Q  |  is  less 
that  Qmax  or  jtj  is  less  than  rmax  .  we  henceforth  require  that  the  weighting 
function  w  is  chosen  so  that 


wn  rn2(!,s)  *  wn  m(l-Q.s-r)  wnim(l,s)  (120) 

for  all  n.m.l.s  and  all  Q  and  r  such  that  j  Q  |  <  Qmax  or  |r|  <  rmax .  That  is, 
the  right  side  of  Equation  120  may  not  change  appreciably  for  values  of  Q  and  r 
for  which  the  amplitude  of  Ry  j  n_m(Q,r)  much  exceeds  zero.  For  this  restriction 
on  the  choice  of  w,  Equations  I  15  and  1 19  can  be  rewritten  3S  follows. 


RT|n,m(|-^.':s-r,s)  =  wnm2(l.s)  RT]nrn(0.t 


(121) 


RT(la“^-la;S3'r,sa)  * 


wn,m'da-rrn,  s3-In)  Rt  |  n.rm’^'T) 


(122) 


a  comparison  of  Equation  121  with  Equation  109  of  Section  3.3  indicates  tnat 
the  non-WSSUS  mutual  coherence  function  for  each  weighted  WSSUS  system 
function  at  any  grid  point  can  be  put  in  time-  and  frequency-  separable  form. 
~he  choice  of  w.  within  the  restrictions  stated  above,  allows  complete  control 
ever  :ne  separable  factors  Jfs)  and  5(1),  and  their  associated  Fourier 
transforms,  and  r^-rp,  discussed  m  Section  5.  These  'our  'actors  are 

nencefortn  denoted  2n(s),  tfm(l),  fn(^-v)  and  fm(£-r()  to  conform  to  the 
notation  of  this  section,  and  to  emphasize  that  one  is  'ree  to  vary  these 
separable  factors  for  each  point  in  the  grid.  Eacn  wn  m(l,s)  snouid  oe  chosen 
equal  to  the  square  root  of  the  desired  product  of  ?n(s)  and  t^O)  at  the  fn,m)in 
grid  point. 


4.2  CHOICE  OF  WEIGHTING  FUNCTION  W  LEADING  TO  PIECE-WISE-LINEAR 
VARIATION  OF  R T  IN  TIME  AND  FREQUENCY. 

R Tda.la;ta-Sa^  wl11  vary  in  a  P'ece-wise-linear  fashion  in  time  and  carrier 
frequency  if  the  time  and  frequency  dependent  factors  of  the  separable  form  of 
pT|n,m  at  sacn  9nd  P°int  takes  the  following  form. 


*n(s) 


_A_(  S/ATTn}  Rect{s/ATTn-0.5)  + 
_A_(  s/ATln)  Rect{s/AT;n*0.5} 


(123) 


3m(D 


where 


_A_{  l/AFTm}  Rect{l/AFTm-0.5}  + 
^V_{  t/AF^}  Rect{l/AFAm*0.5} 


(124) 


ATTn 

= 

Tn+ 1  “  Tn  = 

AF*n*  i 

(125) 

AT^ 

= 

Tn  “  Tn-1  = 

atV, 

(126) 

AFrm 

= 

Fm  + 1  “  Fm  = 

AF^n+ 1 

(127) 

E 

— ► 

U_ 

<] 

= 

"  Fm-l  = 

AFfn-i 

(128) 

ana  _A_(-}  and  Rect{-}  are  defined  below. 

i  -  I  x i  for 

-A_(x)  =  \  11 


for  -1  <  x  £ 
otherwise 


(129) 


Rect(x) 


for  -1/2  <  x  <  1/2 


otherwise 


rne  weighting  function  wnim  ,  resulting  from  the  above 
is  given  by  the  following  equations. 


wn>m(l.s)  =  [  Sn(s)  3m(l)  ]  1/2  =  (i. 

-[  1/2  r  1 

s/ATTn)  Rect{s/ATTn-0.5)  _A_{  l/AFTm)  Rect{l/AFrm-0.5) 

*  X  ♦ 

_A_{  s/ATAn)  Rect{s/ATV0-5}  _A_(  l/AFAm}  Rect{l/AFAm*0.5) 

t  J  l 


1/2  = 


030) 

form  of  ?n(s)  and  ^(1). 

(13D 
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